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The implications of a QCD phase transition at high temperatures and densities for 
core-collapse supernovae are discussed. For a strong first order phase transition to quark 
matter, various scenarios have been put forward in the literature. Here, detailed numerical 
simulations including neutrino transport are presented, where it is found that a second shock 
wave due to the QCD phase transition emerges shortly after bounce. It is demonstrated that 
such a supernova banging twice results in a second peak in the antineutrino spectrum. This 
second peak is clearly detectable in present neutrino detectors for a galactic supernova. 

§1. Introduction 

Core-collapse supernova are enigmatic astrophysical cataclysmic events. The 
observation of supernova SN1987A and the measurements of its emitted neutrinos 
inspired to put forward unconventional explosion mechanisms. The detection of 
neutrinos was puzzling as there was a discrepancy in the timing of the events from 
various neutrino detectors. However, De Rujula pointed out that the measurements 
are consistent with the emission of two separate neutrino bursts delayed by about 
five hours P He concluded that a supernova may bang twice and suggested that a 
collapse to a black hole triggered by accretion can explain the second neutrino peak. 
Hatsuda argued that a hypothetical strange star, a quark star bound by interactions 
only, could have formed within a secondP He pointed out that the released energy 
from the transition to a strange star is comparable to the energy emitted in super- 
novae. The phase transition to quark matter was studied for an adiabatic collapse by 
Takahara and SatoP 1 They argued that the released latent heat increases the core 
temperature drastically thereby generating a prolonged emission of neutrinos. In a 
full hydro simulation with a phase transition a second shock wave emergeo® but neu- 
trinos were not considered in the simulation. Drago and Tambini introduced another 
intriguing scenario where the formation of strange quark matter leads to a prompt 
bounce.^ The mixed phase forms during the supernova collapse, causes a soften- 
ing of the equation of state. At higher densities the equation of state stiffens again 
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for being compatible with the observed pulsar masses. Nakazato, Sumiyoshi, and 
Yamada performed a full supernova simulation including neutrinos and the quark- 
hadron mixed phase.® ^ They started with a rather high initial mass of IOOMq 
finding no second shock wave but only that the softening of the equation of state 
just shortened the timescale for the collapse to a black hole. Here, we report on the 
implementation of the QCD phase transition for core-collapse supernova simulations 
with neutrinos included for lower initial progenitor masses.®"^ A second shock 
wave is generated by the appearance of the quark-hadron mixed phase leading to 
an explosion. When the heated shock material is running over the neutrino sphere, 
a second burst of neutrinos is released in antineutrinos. This second peak can be 
observed with the present neutrino detectors Super-K and IceCube for a galactic 
supernova as shown by Dasgupta et alP^* 

§2. QCD Phase Transition in Neutron Stars 

The QCD phase diagram is largely unknown at high densities and temperatures. 
The asymptotic freedom of QCD ensures that at large enough energy scales matter 
should be described in terms of free quarks and gluons and not by hadronic degrees 
of freedom. At high densities, quark matter could be in a chirally restored but not 
deconfined phase which is dubbed the quarkyonic phased Compact stars consisting 
of pure quark matter, so called strange stars, could be bound just by interactions 
with quite exotic properties as shown within the MIT bag modeP^~^ (see also 
Ref. [To])), In perturbative QCD calculations to 0{offi^^^ the properties of these 
strange stars were found to be surprisingly similar to the results of the MIT bag 
model. Strange stars have quite similar maximum masses and radii compared to 
ordinary neutron stars. However, a strange star is hypothetical as one has to assume 
that strange quark matter is more stable than ordinary nuclear matter. 

Hence, one has to match to hadronic matter at low densities. The onset of the 
mixed phase could be as low as (1 — 2)uq for a large range of values for the MIT 
bag constants, see e.g. Ref. I22|) . Sufficiently high densities are reached in the core 
for a 1.3M Q neutron star to have quark matter present. For a strong first order 
phase transition a third family of compact stars appears in the mass-radius diagram 
besides white dwarfs and neutron starsP^'^ 1 Signals for such a phase transition and 
the appearance of a third family have been discussed intensively in the literature, as 
e.g. the spontaneous spin-up of pulsars^ 1 and an exotic mass-radius relation with 
the so called rising twins!^ 1 Also the collapse of a neutron star to the third family 
can release gravitational waves, 7-rays, and neutrinos. The signal of the QCD phase 
transition for core-collapse supernovae will be the subject of the next section. 

§3. QCD phase transition in core-collapse supernovae 

The conditions for phase equilibrium in supernova matter have been devised in 
detail in Ref. 25). Several different cases have been considered depending on the 
locally and globally conserved charges in the thermodynamic system at hand which 
for supernovae matter would be the proton or lepton fraction in addition to the 
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Coulomb charge. An interesting case concerns the one, where the charge is assumed 
to be locally conserved not globally. Then chargeless bubbles appear in leptonized 
matter at the early stage of the proto-neutron star's life which form a mixed phase 
due to a nonvanishing locally conserved lepton number. The mixed phase disap- 
pears when the neutrinos leave the star which can cause a delayed collapse with a 
pronounced emission of neutrinos Proto-neutron star evolution with quarks have 
been simulated in the work of Pons et alP^ They find that the onset of the quark 
phase in core-collapse supernovae occurs during the late time evolution of the proto- 
neutron star. The timescale for quark matter to appear was found to be typically 
(5 — 20)s which is well after the bounce. The late onset of the quark phase is due to 
a large bag constant, I? 1 / 4 > 180 MeV which was chosen to get a large neutron star 
mass. 

There are several arguments for quark matter to appear at rather low densities 
in astrophysical systems. First in /3-equilibrium, strange quark matter is favored 
in comparison to ordinary nuclear matter due to the additional strangeness degree 
of freedom. Secondly, low values of the proton fraction are also an advantage for 
quark matter due to the large asymmetry energy of nuclear matter (the situation 
gets more subtle when effects of color-superconductivity are considered, see Ref. |28|) 
for more details). Thirdly, the phase transition line is located at lower densities 
at higher temperatures. The total effect is that the production of quark matter in 
supernovae material can occur at quite low densities so that its production at bounce 
is possible. For a temperature of T = 20 MeV the phase transition to the mixed 
phase can be located just slightly above normal nuclear matter density for Y p = 0.3. 
It is important to note that for matter in heavy-ion collisions the phase boundary is 
shifted to much larger values. It is also crucial to realize that the phase transition is 
not necessarily related to deconfinement, so the naive picture of overlapping hadrons 
making the transition can be misleading. In principle, any strong first-order phase 
transition can be envisioned in dense matter for our purpose here, as one related to 
chiral symmetry restoration. 

There are only two supernova equations of state commonly used for the nuclear 
phase, the one of Lattimer and Swesty^ 1 and the one of Shen et alP^ Advanced 
modern equations of state are presently being developed. Medium effects for nuclei 
have been studied in Ref. IHTj) and applied for supernova conditions by Typel et alF^ 
The statistical approach of Botvina and Mishustin for supernova matter utilizes 
methods from nuclear multifragmentation. 33 ' A complete supernova equation of 
state has been developed with an excluded volume which takes into account the 
whole set of nuclei of the nuclear chart in Ref. I34p . In the following we will use the 
equation of state from Shen et al. where the phase transition is added by using the 
MIT bag model. 

It turned out that a particular critical point to check the equation of state is the 
maximum mass for cold neutron stars. Recent mass measurements of the pulsar PSR 
1903+0327 indicate a quite high mass of M = (1.67 ± O.Ol)M PJffl A detailed 
discussion on mass and radius constraints for cold neutron stars can be found in 
Refs. I37|) - f39]) . The mass-radius curve for the equation of state used in the supernova 
simulation is depicted in Fig. Q] for different values of the MIT bag constant. A 
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Fig. 1. The mass-radius relation for cold hybrid stars for different values of the MIT bag constant 
and a s corrections. Present constraints from pulsar mass measurements are for PSR1913+16 
1.44M© and for PSR 1903+0327 l.&7M c ^ m (figure taken from Ref.©). 

lower value results in a higher maximum mass. Here we find maximum masses of 
M max = 1.56M for B 1 / 4 = 162 MeV and Af max = 1.5M for B 1 / 4 = 165 MeV. If 
one includes corrections from one-gluon exchange, the maximum mass increases to 
-^max = 1.67M0 (for a s = 0.3) while the critical density for the onset of the mixed 
phase is tuned to be similar by lowering the bag constant to 1? 1//4 = 155 MeV.®'^ 
The strong first order QCD phase transition has a significant impact on the 
dynamical evolution of core-collapse supernovae. Full hydrodynamical simulations 
including neutrinos have been performed for progenitor masses of M = 1OM and 
M = 15M in Ref. 8). A few hundred milliseconds after the bounce, a second shock 
wave develops due to the formation of the new high-density phase. The mixed phase 
has a lower adiabatic index, collapses and hits the high-density core which has a 
much larger adiabatic index. A shock front forms, travels outwards and accelerates 
at the steep density gradient on the surface of the proto- neutron stars. When the 
heated material passes the neutrinosphere, antineutrinos are released in a burst. The 
antineutrino spectrum will be markedly different from the conventional picture as 
a pronounced second peak appears shortly after the time of collapse. In the sim- 
ulations runs the quark core appears between t p b = 200 to 500 ms depending on 
the progenitor mass and the chosen bag constant. The transition to strange quark 
matter involves the production of strangeness during the conversion process. In 
hadronic matter, hyperons can appear to some extent in supernova material due to 
/3-equilibrium and thermal production, see Refs. I4"0"|). l4ip who find a hyperon fraction 
at bounce (assuming T ~ 20 MeV) of about 0.1%. This small amount of hyperons 
already present enables the nucleation of strange quark matter via fluctuations of 
strangeness.^'® The results are highly sensitive to the surface tension between 
the old and the new phase, which is a basically unknown quantity. The nucleation 
timescales computed are below the timescale of a supernova for critical surface ten- 
sions of a < 20 MeV fm~ 2 , which are usually considered to be small but which are 
not unreasonable. 

The antineutrino peak from the QCD phase transition can be detected with the 
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Fig. 2. The expected signal of antineutrinos for a supernova banging twice at a galactic distance 
of 10 kpc as seen in the Super-Kamiokande detector (left plot) and in IceCube (right plot). A 
spectacular pronounced second peak is clearly seen in the simulated data, here at a time of 
around 260 ms (figures taken from Dasgupta et alP^). 



present neutrino detectors Super-K and IceCube, see Dasgupta et alP3 These detec- 
tors are mostly sensitive to antineutrinos by inverse /3-decay reactions {v e p — > ne + ). 
Adopting the antineutrino spectrum from the supernova simulation, the observed 
antineutrino signal from a supernova would clearly reveal a second peak. Fig. 
shows the expected signal for Super-K and IceCube assuming a galactic supernova 
at a distance of 10 kpc. The second peak in the antineutrino spectra sticks out quite 
strikingly showing the high sensitivity of the present neutrino detectors. Unfortu- 
nately, the sensitivity of the neutrino detectors at the time of SN1987A were orders 
of magnitude lower, otherwise one might have already seen that a supernova can 
indeed bang twice. 
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